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Except for a few special states, computing quantum discord remains a complicated optimization 
process. In this paper we present analytical solutions for computing quantum discord of the most 
general class of A-states and the criteria for each analytical solution to be valid. We discuss param¬ 
eter regions that correspond to different analytical solutions and explain the underlying reasons for 
such structure to exist. We apply our formalism to study both arbitrary A'-states and A-states with 
certain symmetries. We find that our analytical formalism is in excellent agreement with numerical 
calculation of quantum discord in both cases. 

PACS numbers: 03.67.-a,42.50.Ex 


I. INTRODUCTION 


Quantum discord 0,1 as a measure of the quantum¬ 
ness of correlations has lead to revaluation of the ad¬ 
vantage of quantum systems which are affected by de¬ 
coherence. Numerous are the quantum protocols which 
can be performed with quantum discord M instead of 
entanglement Q. As a measure of non-classical corre¬ 
lation, quantum discord is also recognized as a funda¬ 
mental property of quantum systems fTot - TTrl . Although 
the definition of quantum discord is conceptually clear, 
computing it has been difficult since it requires the op¬ 
timization over all possible measurements on one qubit. 
The most general class of quantum states for which quan¬ 
tum discord can be studied anal ytic ally has been shown 
to be the A states 05). A-states[H| are encountered fre¬ 
quently in different contexts of physics 07l - [23| . Although 
there have been a few methods for finding analytical re¬ 
sults for two-qubit A-states (24l - [30| they are not always 
correct. Extensive numerical studies have shown that the 
two known analytical formulas are not always valid and 
errors exist in the calculation of quantum discord with 
such formulas [31 ] . 

In this paper we consider the most general form of 
A-states and we find that the quantum discord for the 
majority of the states can be calculated with two closed 
form formulas. We also determine analytically the cri¬ 
teria for these formulas to be valid and the solution for 
the rare cases where the two formulas are not valid. This 
provides a way to calculate the quantum discord for any 
A-states without performing the optimization process. 
The article is arranged as follows: In Sec. |TT] we intro¬ 
duce the definition of quantum discord 0]. In Sec. m 
we discuss the derivation of our formalism. In Sec. m 
we discuss examples and applications to various A-states 
and the comparison to numerical calculation. Finally, in 
the last section we summarise our principal results and 
future directions. 


II. QUANTUM DISCORD 

Quantum discord arises from the difference between 
two classically equivalent definitions for mutual informa¬ 
tion that are related through Bayes’ theorem 0. Firstly, 
the quantum mutual information, I(pab), is the total 
shared information between two quantum subsystems A 
and B being in a state pab, specifically [32) 

I(pab) = S(pa) + S(pb) — S(pab), (1) 

where S(pa ), S(ps), and S(pab) are the von Neumman 
entropies |3J] of the reduced states pa = Trs(pAB), 
Pb = Tva{pab), and the total state pab respectively. 
Secondly, the quantum mutual information, J (B\A), can 
be defined based on the quantum equivalence of Bayes’ 
theorem. It represents the information gained about the 
subsystem B after measuring the set n a on subsystem A 
and it is given by 

J ( B\A ) = S(p B ) - min S (pb\a) ■ (2) 

Here S (pb\a) = Y, a Pa S (pb\o) is the average condi¬ 
tional entropy of system B , where p a = Tr(Jl a pABA\) 
is the probability to obtain the outcome a and p B \ a = 
IlaPAsK/Pa is the corresponding reduced density for 
subsystem B. The minimization of conditional entropy 
is performed over all the possible positive operator valued 
measure (P.O.V.M.) 0ij. Quantum discord is defined as 
the discrepancy between I [pab) and J (B\ A), 

D{B\A) = I(p AB ) — J(B\A), 

= mia S ( p B \A ) - S(p AB ) + S(p A )- (3) 

Quantum discord can be understood as a measurement 
of the non-Bayesianism of the state pab, or the amount 
of shared information that is not accessible by local mea¬ 
surements. It should be noticed that the information 
about B given a measurement on A is not necessary the 
same as the information of A given a measurement on B 
and therefore D(B\A) ^ D(A\B). The definition given 
in Eq.([3]) shows the challenge of calculating discord is to 
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find the optimal measurement that minimizes the average 
conditional entropy S (pb\a)- Determining the optimal 


To simplify the analysis, we define following functions 


measurement is the focus of our discussion in the next 

B 2 = w 2 + z 2 + 2wzcos(24>), 

(11) 

section. 

A 3 ±Aicos0 

G±= 2 

(12) 

III. FORMULISM 

R± = y]G 2 ± + B 2 sin 2 9, 

(13) 


and constants 


Let us consider two qubits A and B in A-state px, 
which in the logic basis { 00), 01), 10), 111)} has the 

Ai = a — b — c + d, 

(14) 

following matrix representation 

A 2 = a + b — c — d, 



Px = 


( a 0 0 w 

0 b z 0 

0 z c 0 

V w 0 0 d 


( 4 ) 


Here the diagonal elements are complementary probabil¬ 
ities, a + b + c + d= 1, and the positivity of the density 
matrix requires |u;| < yfad and \z\ < y/bc. We consider 
w and z real and positive numbers since their phases can 
be removed by local unitary operation (2] and do not 
change the quantum correlations between A and B. For 
arbitrary two qubit states of rank 2, the P.O.V.M that 
minimizes the conditional entropy is the projective von 
Neumann measurement [23] ■ For states of rank 3 and 4 
F. Galve et al. [291] showed that projective measurements 
are almost sufficient, except for states that appear with 
probability 10 -2 and in which case the deviations are on 
average of the order of 1CH 6 . The projective measure¬ 
ment bases are characterized by the angles 6 and (j) of 
the Bloch sphere, as follows 

|+) = cos(0/2)|O) + e*^sin(0/2)|l), (5) 

|—} = sin(0/2)|O) — e l(t> cos(0/2)|l), (6) 

with 9 G [0,7r] and <f> G [0, 2n\. We consider measuring A 
with the above measurement bases, n± = |±)(±|. The 
corresponding outcome probabilities p± are 


P± = t ;(1 ± A 2 cos9), 


( 7 ) 


and the average conditional entropy of B as a function 
of 6 and q.5, F(9 , <j>), is written as 


F{9,(/)) = p + S(p B \+) + p_S(p B \-), 

= C^+ + 7 1 - + Q+ + Q -), 


with 


T± {9 1 (f>) = (jp+ ± R+) log 2 
Q± (0 ,<t>) = (P- ± 7?-)log 2 


P+ ± 7^4 
2 P+ 

P- ± R- 
2 p- 


( 8 ) 


( 9 ) 

( 10 ) 


A$ = a — b + c — d. 

It is worth noting that F(9,(j>) = F(9,2n — cj>) and 
F(9, 4>) = F(tt — 9, </>), i.e. F is symmetric with respect to 
(f> = 7r and 9 = tt/2. It is therefore sufficient to consider 
4> € [0,7r[ and 9 G [0, 7r/2]. Furthermore, the only (j> de¬ 
pendence of F is in the term wz cos(2<(i) of B 2 in Eq. lfTTl) 
and F becomes independent of (f> when w or z is zero. 
We also note Ai, A 2 and A 3 are the expectation values 
of the ^-component of the Pauli vector, A\ = (o-^crf), 
A 2 = (erf) and A 3 = (af )[33j. The minimum of F(9,(F) 
can be found by using the extreme value theorem 1351 . 
The critical points c = ( 9 C , <j> c ) are determined by solv¬ 
ing the first partial derivatives equations dF/89 = 0 and 
8F/d(f> = 0 and whether the critical point c corresponds 
to a minimum is determined by the second partial deriva¬ 
tives test. i.e. the second partial derivatives and the 
Hessian matrix determinant have to be simultaneously 
positive at c 


d 2 F 


89 2 


> 0 , 


d 2 F 


dej) 2 


> 0, det 


d 2 F 

8 2 F 

89 2 

dOdcj) 

d 2 F 

8 2 F 

dOdcf) 

d<p‘ 2 


> 0. 


Otherwise the critical point c can correspond to a 
local maximum, a saddle point, or if Hessian matrix 
determinant is zero, higher order test must be used to 
determine the nature of the critical point. 


We find the first partial derivatives of F as 
8F 


dej) 

8F 

~89 


= —wzsva 2 9su\(2(j))C ( i > , 

sin 9 
~\ ^ 0 ? 


with 

C^ = 

and 


1 p-+R- 

TC l0& TTV-rT 


_L, og 5t±i. 

R+ 8 p+- R+ 


(15) 

(16) 

(17) 


^ AiG- + 2B 2 cos9 p-+R- AiG+ - 2B 2 cos9 , p+ + R+ 

Ce = - 77 -log 2 - - - - -log 2 ■ 


R- 


' p- — R- 


R~ 


p+- R+ 


A 2 log. 


p 2 + (pL - R -) 
P 2 . (p\ - R \)' 


(18) 
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From these expressions we note several key features about 
the critical points: 

• Because 0 < R± < p±, C^ is always positive and 
the partial derivative of F with regard to </> in 
Eg. (1151) is zero only when sin(2<^>) is zero for any 
9, i.e (f> c = 0, 7 t/ 2 , and when 0 = 0 for any (j). We 
note that if 9 is zero then the minimization does 
not depend on (f>. At these <j) c the second derivative 
d 2 F/d<j) 2 becomes 


d 2 F 

d<j) 2 


= 2wzC ( f >= o sin 2 9 > 0 

( 0 , 0 ) 


V 9 7 ^ 0 and 


d 2 F 

d(f> 2 


(«■*) 


—2 wzC<t,= * sin 2 9 < 0 


(19) 


( 20 ) 


V0^O. 

Since for any 9 the second derivative d 2 F/dcf> 2 is 
always negative when </> = 7r/2 the minimization 
problem becomes one variable minimization prob¬ 
lem where (j> c = 0. As F(6, 0) becomes a one vari¬ 
able function is not necessary consider the Hessian 

I 


determinant it is sufficient study the behavour of 


• The derivative Eq. (usd is zero at either sin 9 = 0 or 
Cg = 0. The former equation provides 9 C = 0. The 
equation Cg = 0 has the obvious root 9 C = 7 t /2 and 
a special root 9 e . This special root depends on the 
density matrix elements and only appears in some 
very particular states, such as shown in FigOJh). 
The angle 9 e is not considered in previous analytical 
studies [III HH, |26[. It corresponds to the “error 
cases” reported in numerical calculations [31]. 

Based on these features, the optimization problem re¬ 
duces to studying the second derivative of F{9, 0) with 
respect to 9 evaluated at the critical angles 9 C = 0, 7 t/ 2 , 
and 9 e . We notice that the measurement bases become 
a z for 9 C = 0 and <r x for 9 C = 7 r/ 2 . The value for 9 e 
remains dependent on density matrix elements and we 
denote this measurement as er e - 

The second derivative evaluated in the three 9 C de¬ 
pends on the behaviour of two quantities Co and C + given 

by 


C 0 = A 2 log 2 


cd (1 + A 2 ) 2 
ab (1 - A 2 f 


. , dd _ , x 2 

A 1 log 2 -—I- 2 (w + z) 
be 



, a 
lo §2 T + 


1 

c — d 



and 

C+ = 4 r (A 1 A 3 — 4 A 2 r 2 )“ — 4 (w + z) 2 (l — 4?' 2 ) (4r 2 — A\) In — 


( 21 ) 


( 22 ) 


where r = A 2 + 4 (w + z) 2 /2. The behaviour of Cq and 
C + determines which of F(0,0), F(n/2,0), and F(6 e , 0) 
is a minimum. The relationship is summarized in Table 
|T| It is worth noting that for more than 99% parame¬ 
ter space [§l|, F(0,4 ») is minimized either at F( 0 , 0 ) or 
F(7r/2,0), which can be calculated as follows, 


F{ 0 , 0 ) = —alog 2 ——61og 2 —— clog 2 -^—-dlog 2 

a+b a+b c+d c+d 


and 


(23) 


’(?•»)■ 


1+\/A^+4(w+z) 2 l+\/Al+4(w+z) 2 

o ^§2 o 


1-\/A 2 +4{w + z) 2 1-y/Ag +4(«’ + a) 2 ( 2 1 ) 


The analytical solution also agrees with the particular 
expressions obtained in Ref.|26[ for symmetric X-state, 
i.e with b = c in the density matrix of Ecp(U). We find 
that there are three possible minima for an arbitrary X- 
state and we find for the first time the conditions for 


special cases a e and <7?. In the following section, we dis¬ 
cuss examples where different measurements are neces¬ 
sary and use our formula to calculate the quantum dis¬ 
cord. 


IV. RESULTS 

A. General X-states 

First we consider two examples of general X-states 
where the rare measurement cases exist. In Figure HJa), 
we show an example where the optimal measurements 
changes from <j z , cr? and <r x as w and z changes with the 
diagonal density matrix elements fixed. For the parame¬ 
ter space where w and z go from zero to approximately 
0.0502681 the optimal measure is always oy, as indicated 
by the dashed line. For the parameter space where w and 
z go from approximately 0.0540158 and above, the opti¬ 
mal measure is a x , as indicated by the solid line in the 
plot. In the extremely small region (shown as the insert 
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TABLE I. Criterion for different optimal measurement (M) in 
terms of Pauli matrices and the solution to min F(9, <j>). The 
measurements of a z and a x are the predominant cases for an 
arbitrary X-state. They include more than 99% of the pa¬ 
rameter space. For the case where the measurement is listed 
as any, F(9,(j>) is the same for any projective measurement. 
For the case where the measurement is listed as <r e , it is nec¬ 
essary find 9 e from Cg = 0 solving for every set of density 
matrix elements. For case of a?, min F{9,(j>) is the minimum 
of F(0,0) and F(f,0). 


M 

Condition 

min F (9, (f>) 

any 

Co 

= 

0VC+ 

= 

0 

F( 0,0) 

(Tz 

Co 

< 

< 

o 

> 

0 

F{0, 0) 

O' x 

Co 

> 

o 

> 

■P 

< 

0 

0) 

O e 

Co 

> 

o 

> 

p 

> 

0 

F(9e, 0) 

O'? 

Co 

< 

0AC+ 

< 

0 

min \ F(0, 0), F(^, 0) j- 


plot), where w and z go from approximately 0.0502681 
to 0.0540158, the measurement is denoted as cr?. This is 
the region where both Co and C+ are negative and the 
optimal measurement has to be determined by calculat¬ 
ing F(0,0) and F(ir/2,0). In Figure QJ&), we show an 
example where the optimal measurements changes from 
a z , a e and a x as w and z changes with the diagonal 
density matrix elements fixed. For w from zero to ap¬ 
proximately 0.233563 and z from zero to approximately 
0.111803 the optimal measure is always a z , as indicated 
by the solid line. When w is greater than approximately 
0.116136 and z is greater than approximately 0.111803, 
the optimal measure is a x , as indicated b the dashed 
line in the plot. In the small region between the dashed 
and solid lines (shown in the insert plot), neither F( 0, 0) 
nor F(tt/2,0) is the minimum. In this region the op¬ 
timal measurement is given by the special measurement 
<r e . We also note that our criterion reproduces the results 
obtained in references m and [28|. 

Figure (j2j shows a comparison between numerical [.36j| 
and analytical results for quantum discord. Here dis¬ 
cord is calculated for the density matrices discussed in 
Fig©6) with w = 0.05, and the optimal measurement 
can be either a z , a x or a e . The numerical solution of 
discord is given by the dashed line. When z goes from 
zero to z 0 ~ 0.0655096 the optimal measure is a z (red 
line). When z goes from z + ~ 0.0661362 to its maxi¬ 
mum value Zmax = Vbc the optimal measure is a x (green 
line). When z goes from zo to z+ the optimal measure 
is <r e (blue dots). In this last case it is necessary solve 
Ce(6 e , 0) = 0 for every z. The numerical and analytical 
results agree. In Fig.©) we show 6 e as a function of z in 
the special measurement regime as indicated in Fig. (f2|) . 
Although in Fig.© the transition between the optimal 
measures a z and a x appears sharp, Fig. d3j) shows that for 
z from zo to z + , the special angle 0 e goes continuously 
from 6 = 0 to 6 = 




FIG. 1. Regions where the different cases arise as function 
of w and z, for two particular X-states: (a) the diagonal 
elements are a = 0.5, b = 0.3, c = 0.1 and d = 0.1, ( b ) the 
diagonal elements are a = 0.0783, b = 0.1250, c = 0.1 and 
d = 0.6967. The solid and dashed line correspond to Co = 0 
and C+ = 0, respectively. 



FIG. 2. Quantum discord, D(B\A), as function of z for 
a = 0.0783, b = 0.1250, c = 0.1, d = 0.6967 and w = 0.05. 
Black dashed line correspond to numerical solution of Eq. © 
solved for every value z. Red and green lines are our analyt¬ 
ical calculation based on Eg. (12311 and Eg. (1241) . The optimal 
measurement for the red line is a z . The optimal measurement 
for the green line is a x . The blue dots are the result from the 
special measurement <j e , which is determined by the roots of 
Cg = 0 for a given z value. 
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FIG. 3. Special angle 9 e of the measurement a e as a function 
of z in the special measurement region of the state given in 
Fig. 2. In this region, the special changes continuously from 
zero to 7 t/2. 


B. States with symmetries 



FIG. 4. Quantum discord for two spins in the XXZ model 
with a = 0.1 and w = 0. Black dashed line correspond to 
numerical calculation. The red line corresponds to where the 
optimal measurement is u z and the quantum discord is de¬ 
termined through calculating Ea. (l26l) . The green line cor¬ 
responds to where the optimal measurement is a x and the 
quantum discord is determined through calculating Eq. m- 


The analytical derivation in Sec. also shows that the 
behaviour of F strongly depends on the constants Ai as 
defined on Eq. (fl5l) . These constants corresponds to the 
expectation values of a z s and they become zero when 
certain symmetries are present. For example, a system 
with the spin-flip symmetry, such as the XXZ model, 
satisfies A 2 = {af) = 0 , A 3 = (erf) = 0 and A\ = 
(<J^u z ) = 4a — 1 with 0 < a < 1/2. The most general 
form of the reduced density matrix of the two spins in 
such systems are written as 


f a 0 0 w \ 

0 5 — a z 0 

0 z \-a 0 ’ 

yw 0 0 a ) 


(25) 


with 0 < w + z < 1/2. For such systems the only optimal 
measurements are <j x and a z . For 0 < + z < 12a — ^ |, 

the optimal measurement is <r x and the minimum of 
F{0, (f>) is 


F( 0 , 0 ) = - 2 a log 2 ( 2 a) - (1 - 2 a) log 2 (l - 2 a). (26) 

For 12a — i| < w + z < i, the optimal measurement is 
a z and the minimum of F is 

F(/o)=- 1+2t ; +,| i, e , 1+2( ; +;| 

— 1 ~ ~ ~ 1°S2 1 ~ + ~ ■ (271 

As a example we consider the state of two spins in the 
XXZ model, which corresponds to w = 0 in Ea. (l?5l) . In 
Fig-dl]) we show the typical behaviour of the quantum 
discord between two spins in the XXZ. We consider 
varying z and a = 0.1 in the density matrix. The opti¬ 
mal measurement changes from cr z to a x as z increases. 



FIG. 5. Optimal measurements on the z — a plane with w — 0. 
This corresponds to the XXZ model. Positivity of the density 
matrix requires that z < 1/2 — a. The light grey area denotes 
the parameter region for the optimal measurement a z and 
the dark grey area for the optimal measurement a x . The 
dashed line corresponds to z = \2a — 1/2 1 and the solid line 
correspond to z = 1 /2. On the dashed and solid lines Co and 
C+ are both zero. This corresponds to the first case in table 
0. where F is independent of the measurement angle. 


On the transition point from tr*. to cr x , the discord is 
independent of the measurement angle. We find that nu¬ 
merical and analytical calculations agree perfectly for all 
z. In Fig. m we show the parameter regions for optimal 
measurements a z and a x for a general two-spin state in 
the XXZ model. We find that when z is greater than 
1/2 — 2 a and less than 2 a — 1 / 2 , the the optimal mea¬ 
surement is a x . In the rest of the parameter space the 
optimal measurement is a z . On the border of the param¬ 
eter regions, as indicated by the dashed and solid lines, 
quantum discord does not depend on the measurement 
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angle. 

V. SUMMARY 

In this paper we analyse the quantum discord of arbi¬ 
trary two-qubit X-states. We describe a straightforward 
criterion for determining analytical solution to quantum 
discord and it has shown excellent agreement with nu¬ 
merical calculations for both arbitrary X-states and X- 
states with special symmetries. We find that the op¬ 
timization of measurements to minimize the averaged 
conditional entropy, the key and most computationally 
intensive element of calculating quantum discord, is re¬ 
duced to three possible cases. Two of them are the op¬ 
timal measurements discussed in previous literature and 
one additional one is the ’’error” cases that have been 
observed numerically previously. We provide a complete 
framework to categorize these measurements and deter¬ 
mine their parameter regions based on the analytical cri¬ 


teria we derived. In addition to simplifying the quantum 
discord calculation, our analytical expressions reveals the 
relationship between the optimal measurements and sym¬ 
metries of the states. 

Our analytical expression also has experimental im¬ 
plications in a sense that it provides a way to calculate 
the quantum discord based on the expectation value of 
spin components, which can be measured experimentally. 
Whether our analytical solutions provides more insight 
on the relationship between the behaviour of the opti¬ 
mal measurement and the critical phenomena, as previ¬ 
ous numerical study suggests [? ], remains to be further 
investigated. 
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